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Abstract. 

Hadronic emission from parsec size radio lobes in active galactic nuclei (AGN) is discussed. The lobes are composed 
of shocked jet plasma and expected to be filled with high energy particles. By using the Monte Carlo simulation, 
we calculate the photon spectra from the lobes including photo-meson interaction processes. When the synchrotron 
emission from primary electrons is bright, synchrotron-self-Compton component is dominant in 7-ray bands. The 
hadronic emission from the lobes can be dominated in 7-ray bands when the primary emission is not very bright. 
Proton synchrotron component arises at sub MeV band. The synchrotron emission radiated from secondary pairs 
produced via photo-meson cascade emerges in ^GeV-TeV energy ranges. These high energy emission signatures 
provide a test for proton accelerations in young AGN jets. 
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1. Introduction 

The progress of VLBI (Very Long Baseline Interferometry) 
observations reveal the existence of compact radio- loud ac- 
tive galactic nuclei (AGN) with its linear size LS < 1 kpc 
(e.g., Readhead et al. 1996; O'Dea & Baum 1997). These 
compact radio sources are currently understood as young 
progenitor of large radio galaxies (e.g., Fanti 2009 for re- 
view). Furthermore, recent investigations discover smaller 
radio lobes with LS - 0(10 pc) (e.g., Snellen et al. 2004; 
Orienti et al. 2008). It is well known that recurrent radio 
sources also possess mini lobes inside them (e.g.. Walker 
et al. 2000). High energy emissions from these compact ra- 
dio sources have been recently explored by some authors 
(e.g., Stawarz et al. 2008; Kino et al. 2009). However, the 
previous works focus on the leptonic emissions and little 
is known about a hadronic emission in them. 

Since mini radio lobes with LS ^ 10 pc are closely 
located near the AGN core, the lobes are expected to be 
in dense external radiation fields (Fig. 1). Furthermore, 
together with electrons, the shocks accelerate protons as 
well. Therefore high energy protons are naturally expected 
to be filled with the radio lobes although its total amount 
is still unknown. Therefore p^y interaction is inevitable for 
the mini lobes. Here, we examine interactions between 
high energy protons and surrounding target photons. As 
for external photons, we focus on UV photons from the 
standard accretion disks. 



2. Model 

Following the standard picture of radio lobe formation in 
AGNs (e.g., Begelman et al. 1984), here we briefly review 
the standard picture of radio lobes. When a jet interacts 
with surrounding ambient matter, most of its kinetic en- 
ergy is dissipated via strong shocks. The termination point 



at the tip of the jet is called as hot spots. The hot spot 
is identified as the reverse shocked region of the deceler- 
ated jet. The shocked jet plasma leak from the hot spot 
and they form radio lobes. Therefore we can say that the 
radio lobes are the remnant of the decelerated jets which 
contain relativistic particles. In this work, we examine the 
case when the lobe contains relativistic protons. 

2.L Proton acceleration and coolings 

Hot spots in powerful radio galaxies are one of the 
promising sights for proton acceleration (e.g., Rachen and 
Biermann 1993). Here we show that hot spots in mini lobes 
are also plausible sights for proton acceleration. Hereafter 
the maximum energy of shock accelerated protons is de- 
noted as £^p,max = 1p^p(?' The proton acceleration time 



scale at the hot spots tp^acc,hs 
as 



(,pEp/ [eB\^^c) is estimated 
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where 5hs, and are the magnetic field strength, and the 
ratio of the acceleration timescale and Larmor timescale 
for protons at the hot spot, respectively. Although there is 
no direct estimate for ^hs, the magnetic field strength in 
mini lobes {B) have been recently estimated as 5 ~ 10 — 
100 mG by VLBI observations (Orienti et al. 2008). Since 
the magnetic field strength at the hot spot 5hs should be 
comparable or larger than B. Here we set 5hs ~ 100 mG. 
The value of is a free parameter. For blazars, it is con- 
strained as > 10 (Aharonian et al. 2002). There is little 
constraint on for hot spots in mini lobes and we examine 
the case of = 100 in this work. The timescale of pro- 
ton synchrotron is tp,syn,hs = {QTrrripC^) / {armlEB'^) ^ 
1 X 10^{Ep/10^^ eV)-i(5hs/0.1 G)-^ The high energy 
protons escape from the hot spots via sideways expan- 
sions and they are injected into the mini lobes. When the 
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escape velocity from the hot spot is ^ 0.3 c, the escape 
timescale is accordingly tesc,hs ^ 3 (i?hs/10^^ cm) years. 
Therefore, the typical £^p,max is obtained by the relation 

^p,acc,hs ~ ^esc,hs ^^nd it is 

W-lKlO-(i5;)"'(^)(^)eV(2) 

The Larmor radius of high energy protons with Ep ^ 
10^^ eV is sufficiently smaller than i?hs- Hence, the condi- 
tions of proton acceleration are satisfied. 

Next, high energy protons escape from the hot spot via 
sideways expansions and they injected in the mini lobes. 
Then, they undergo various coolings in the lobes, i.e, adia- 
batic loss, p7 interaction, and proton synchrotron cooling. 
Below, we calculate these processes with the Monte Carlo 
simulation. 

3. Monte Carlo simulation 

To calculate broadband photon spectra, the Monte Carlo 
simulation has been performed in this work. The numeri- 
cal code has been developed by Asano et al. (2008, 2009) 
and references therein. Therefore we do not repeat it here. 
Details are shown in Kino and Asano (2010). Distributions 
of particles and photons are assumed to be isotropic. 
We include the following physical processes: (1) photo- 
pion production from protons and neutrons (p + 7 ^ 
p/n + 7r°/7r+), (2) pion decay (tt^ 27), (tt^ fi^ u) 
and muon decay (/i^ ^ + z/), (3) photon-photon pair 
production (7 + 7 ^ e+ + e~), (4) Bethe-Heitler pair pro- 
duction (j9-h7 ^ p + e+-he~), (5) synchrotron and inverse 
Compton processes of electrons/positrons, protons, pions, 
muons with Klein-Nishina cross section, (6) synchrotron 
self-absorption for electrons/positrons, and (7) adiabatic 
expansion loss. 

4. Results 

Fig. 2 shows the photon spectra for the case of Le = 
1 X 10^^ erg s~^, 1 X 10^^ erg s~^, 1 x 10^^ erg s~^, 
1 X 10^^ erg s~-^, and 1 x 10^^ erg s~^ where Lg is the in- 
jection power of non- thermal electrons. For electron Gyro- 
factor, we set = = 1 X 10^. The injection index is 
assumed as 5 = 2 both for electrons and protons. Based 
on VLBI observations, we set the following parameters 
as = 2 pc, Rhs = 0.3 pc, B = 0.1 G, Vexp = 0.1 c, 
LS = 10 pc, and 7p,min = 7e,min = 10. As for the lumi- 
nosity and temperature of the standard accretion disk, we 
assume I/disk = 3 x 10^^ erg s~-^, and kT^isk = 10 eV (e.g., 
Ostorero et al. 2009). We impose a injection power of high 
energy protons Lp being smaller than total kinetic powers 
of powerful radio galaxies ~ iq^''"^^ erg s~^ (e.g., Ito et 
al. 2008). We assume Lp = b x 10^^ erg s~^ in this work. 

Thin solid lines in Fig. 2 shows the photon spectra 
without proton injection (i.e., Lp = 0). The synchrotron 
spectra from the primary electrons are well studied by 
VLBI observations. The turnover frequency is mainly 
caused by synchrotron self absorption (SSA) (e.g., Snellen 



et al. 2000). According to the synchrotron luminosity of 
the mini lobe Lgyn, the observed turnover frequencies of 
mini lobes show i^ssA ^0.1 — 10 GHz. At the same time, 
primary electrons suffer from fast synchrotron cooling. 
Therefore a break frequency due to the synchrotron cool- 
ing typically appears comparable to or below i^ssA- This 
implies that Lg ~ i^syn- 

Thick solid lines in Fig. 2 represent the photon spec- 
tra with proton injection. It is found that the amount of 
primary electrons controls whether hadronic emission is 
visible or not. For the lobe with Lg = 1 x 10"^^ erg s~^, 
the leptonic emissions overwhelm the hadronic emissions 
at all energy domains. As a result, the spectrum with pro- 
ton injection is almost the same as the one without protons 
injection. Therefore, the bright lobe is not suitable for test- 
ing the existence of high energy protons. As Le decreases 
to Lg < 1 X 10^^ erg s~^, thick and thin lines become sepa- 
rable at 7-ray energy range because synchrotron emission 
from secondary electrons is added in the thick lines. For 
the lobes with Lg < 1 x 10^^ erg s"-*^, the proton syn- 
chrotron bump appears in sub-MeV range with the peak 
at 

'^...„-0.7(^)'(^)MeV. (3) 

In TeV energy range, it is shown that TeV-7 photons are 
significantly attenuated due to soft photons in the lobes. 

5. Summary 

Hadronic emission at extragalactic mini radio lobes is dis- 
cussed in this work. We calculate the broadband pho- 
ton spectra taking interaction and proton synchrotron 
cooling into account. The powers of standard accretion 
disk and injected protons are Ldisk = 3 x 10^^ erg s"-*^, 
Lp = b X 10^^ erg s~-^, respectively. 

For bright lobes with Lg ~ Lgyn ~ 10^^ erg s~^, the 
predicted VHE emission is detectable with Fermi/L AT, 
and HESS. Since it is overwhelmed by the leptonic com- 
ponents in this case, it is hard to figure out whether the 
7-ray emission is hadronic- or leptonic-origin. For the lobe 
with Lg ~ 10^"^ erg s~^, secondary pairs radiate TeV 
7 synchrotron emission and it is detectable by the current 
TeV 7 telescopes. 

For less luminous radio lobes Lg ^ Lgyn < 
10^^ erg s~^, the proton synchrotron bump appears at sub 
MeV energy band. The synchrotron bump from secondary 

pairs created by /i-decay appears at GeV/TeV ranges. 
If the double bumps due to these hadronic processes are 
detected, it would be a support for the hadronic model 
since leptonic model does not easily account for such dou- 
ble bumps. 

6. Future prospects 

We add to comment on the recent radio observation of 
mini lobe 3C 84 which is associated with NGC 1275 
{z = 0.0176). It shows the outburst around 2005 and the 
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new component C3 emerges with the size smaher than 
1 pc (Nagai et al. 2010). Fermi/L AT also detect GeV 7- 
ray emission in the same period (Abdo et al. 2009). Fig. 3 
displays the actual multi-epoch images of 3C 84 obtained 
by VERA observations at 22 GHz (Nagai et al. 2010). 
While the flux of C3 component measured by VERA in- 
creases, GeV 7-ray flux over the period August 2008 - 
August 2009 seems to remain constant (Kataoka et al. 
2010). The difference of the measured flux variations be- 
tween VLBI and Fermi/LAT could indicate a breakdown 
of the one-zone synchrotron-self Compton model proposed 
in Abdo et al. (2009). Hadronic emissions could partly con- 
tribute in 7-ray energy bands, although we did not deal 
with the specific case of 3C 84 in this work. In order to 
nail down the real 7-ray emitter, continuous observations 
by VLBIs are quite important. 

In the near future, the project of VLBI Space 
Observatory Programme-2 (VSOP-2) with the high an- 
gular resolution (Tsuboi et al. 2009) will play a unique 
role for mini lobe observations. Observations will be made 
with the 9-m antenna in orbit together with ground ra- 
dio telescopes, and will achieve the angular resolutions of 
40, 80, and 210 microarcsec at 43, 22, and 8 GHz, respec- 
tively. Especially, observations at 8 GHz will unveil fine 
structures of unresolved mini lobes. Furthermore, the era 
of VSOP-2 will overlap the duration in which the next 
generation TeV 7-ray telescope CTA would be in opera- 
tion (http://www.cta-observatory.org/). CTA will have a 
factor of 5-10 improvement in sensitivity in the current 
energy domain of about 100 GeV to some 10 TeV and an 
extension of the accessible energy range well below 100 
GeV and to above 100 TeV. A future collaboration be- 
tween VSOP-2 and CTA will be one of the best ways to 
explore the origin of 7-ray emission from mini radio lobes. 

Acknowledgements. We are indebted to H. Nagai for providing 
us the image in Fig. 3. 
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Fig. 1. A cartoon of mini radio lobes immersed in ra- 
diation field due to the accretion disk. In this work, 
we examine the case of R = 2 pc, LS = 10 pc and 
Ldisk = 3x10^5 erg s-^. 
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Fig. 2. Predicted broadband spectra from the mini radio 
lobe for Lp = bx 10^^ erg s~^ and Luv = 3 x 10^^ erg s~^. 
Black, green, yellow, red, and blue lines are the case for 
Le = 1x10"^^ erg s'^^lxlO^"^ erg s-\ 1x10"^^ erg s-\ Ix 
]^q42 ^ ^ -|^q41 respectively. The lobes 

are located at the distance D = 100 Mpc. The proton 
synchrotron bump is predicted at ~ sub MeV for smaller 

Le. 
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Fig. 3. Multi-epoch VERA images of 3C 84 at 22 GHz. 

The CI component is assumed to be the core. The C3 
component is the new born one and it proceeds to the 
south. The size is still less than 1 pc. 
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